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We analyze the classical phase space of the hydrogen atom in crossed magnetic and circularly 
polarized microwave fields in the high frequency regime, u sing the Chirikov resonance overlap 
criterion and the renormalization map. These methods are used to compute thresholds to large 
scale chaos and to ionization. The effect of the magnetic field is a strong stabilization of a set of 
' invariant tori which bound the trajectories and prevent stochastic ionization. In order to ionize, 

' larger amplitudes of the microwave field are necessary in the presence of a magnetic field. 
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I. INTRODUCTION 



The chaotic ionization of the hydrogen atom in a variety of external fields is a fundamental problem in nonlinear 
dynamics and atomic physics. Early work, in particular, focused on the ionization mechanism in a linearly polarized 
(LP) microwave field g]. This problem is noteworthy because it showed the general applicability of the "Chirikov 
Resonance Overlap Criterion" ^ to real quantum mechanical systems. This empirical criterion predicts the value 
at which two nearby resonances overlap and large scale stochasticity occurs. There is no doubt that the Chirikov 
criterion is generally robust and, therefore, it is no surprise that Chirikov's criteria has been tried as a way to 
' predict transitions to chaos and ionization dynamics in more complicated circumstances, e.g., for the hydrogen atom 
in circularly polarized (CP) microwave fields. However, the success of the Chirikov criterion in describing the LP 
problem stands in contrast with what seems to have been a somewhat mixed performance when applied to the hydrogen 
atom in rotating microwave fields. Partly for this reason, a good deal of controversy has surrounded the ionization 
mechanism and also the interpretation of the Chirikov criterion when applied to ionization in rotating fields. In this 
article we propose to examine the ionization of the hydrogen atom in a circularly polarized microwave field (CP). 
Our analysis will use an advanced method, the renormalization map |5|, so as to include higher order resonances 
^.!^ ' beyond the Chirikov approach. 

, Interest in the CP microwave problem began with experiments by the Gallagher group ||^ which revealed a strong 

• dependence of the ionization threshold on the degree of polarization. They explained their CP results proposing that 
, in a rotating frame, ionization proceeds in roughly the same way as for a static field, i.e., a static field has the same 

■ effect whether its coordinate system is rotating or not. In a Comment, Nauenberg Q argued that the ionization 
mechanism was substantially more complicated and that the effect of rotation on the ionization threshold had to be 
taken into account. Fully three-dimensional numerical simulations by Griffiths and Farrelly were able to provide 
quantitative agreement with the experimental results. Griffiths and Farrelly and Wintgen 1^ independently 
proposed similar models for ionization based on the Runge-Lenz vector. 

^ , Subsequent theoretical work can be divided into three main categories: (i) classical simulations |10 , (m) quantum 

• simulations, and (Hi) resonance overlap studies. The two main papers on resonance overlap are by Howard who 
rN [ published the first such study of this system, and by Sacha and Zakrzewski In both cases, the Hamiltonian 

■ was written in an appropriate rotating frame, a choice of 'zero-order' actions was made and the Chirikov machinery 
invoked. Somewhat surprisingly, the paper by Sacha and Zakrzewski disagrees in a number of key areas with the 
results of Howard [TT| as well as with some conclusions drawn from numerical studies [p^ . In a paper by Brunello et 
al. , it was shown both numerically and analytically that the actual ionization threshold observed in an experiment 
must take into account the manner in which the field was turned on; essentially, in some experiments electrons are 
switched during the field "turn on" directly into unbound regions of phase space. That is, the underlying resonance 
structure of the Hamiltonian is almost irrelevant because ionization is accomplished by the ramp-up of the field. For 
this reason, the application of resonance overlap criteria in rotating frames can be quite intricate and this provides 
some explanation for the apparent inadequacy of the Chirikov criterion - if ionization occurs during the field turn-on 
time then, of course, resonance overlap is irrelevant. 

In this article, we study the hydrogen atom driven by a CP microwave field together with a magnetic field lying 
perpendicular to the polarization plane (CPxS). The magnetic field has been introduced to prevent ionization in the 
plane. This provides the opportunity to eliminate difficulties associated with the turn on of the field, thus opening 
up the way to a study of resonance overlap in a more controlled manner in the rotating frame. As noted, without the 
added magnetic field all the electrons may have gone before the resonances have had a chance to overlap. 
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The paper is organized as follows: Section II introduces the Hamiltonian and provides a discussion of resonance 
overlap in the rotating frame in both the Chirikov and renormalization approaches. In order to obtain qualitative and 
quantitative results about the dynamics of the hydrogen atom in CPxB fields, we compare numerically, in Sec. Ill, 
Chirikov's resonance overlap criterion with our renormalization group transformation method. The Chirikov method is 
found to provide good qualitative results and is useful because of its simplicity. The renormalization transformation is 
used first to check the qualitative features obtained by the resonance overlap criterion, and to obtain more quantitative 
results about the dynamics by expanding the Chirikov approach. Conclusions are in Sec. IV. 

II. HYDROGEN ATOM IN CPxB FIELDS 

We consider a hydrogen atom perturbed by a microwave field of amplitude F and frequency uj f , circularly polarized 
in the orbital plane, and a magnetic field with cyclotron frequency coc- The Hamiltonian in atomic units and in the 
rotating frame of the microwave field is |l5| : 

H{p,,Py, X, y) = ^^-1^ - J. - n{xpy - yp,) + Fx + ^{x' + y^), (1) 

^ y/x'^ + y^ o 

where VI — Uf — ujc/'i. Following Ref. jl^, we rewrite Hamiltonian (|]) in the action-angle variables (J, L, 9, -0) of the 
problem with ujc = 0. The angles 9 and '0 are conjugate respectively to the principal action J and to the angular 
momentum L. The Hamiltonian becomes : 

H{J,L,9,^P) ^ Ho{J,L,9) + FJ^ J2 t4(e) cos^ + 0), (2) 

n— — oo 

where the integrablc Hamiltonian Ho is 

Ho{J,L,9) = -^-nL + J^^ ^ V„,V^-nC0sn9. (3) 



The coefficients Vn of the expansion are the following ones 
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Voie) = -y, (4) 



Vnie) = - 
n 



Jn{ne) H Jn{ne) 



, for n 7^0, (5) 



where Jn is the nth Bessel function of the first kind and J'^ its derivative. The parameter e is given by 

^2x1/2 



The positions x and y are given by the following formulas : 

+ c>o 



X = ^ K (e) cos{n6 + tp) , 

n— — oo 

+ 00 

J2 Vn{e)sm{n9 + tl;). 

71 — — oo 

The Hamiltonian (^) can be rescaled in order to withdraw the dependence on fl. We rescale time by a factor 57 [ 
we divide Hamiltonian by il]. We rescale the actions J and L by a factor A = Q^^^, i.e. we replace H{J, L,9,^) 
by AiJ( J/A, L/A, 6*, -0). We notice that this rescaling does not modify the eccentricity e. The resulting Hamiltonian 
becomes : 

+ 00 

H = Ha + F'J^ T4cos(n6l + 0), (6) 



n— — oo 
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where 

^°""2P + ^ K^14^-„cosn0. (7) 

n,m— — oo 

and F' — Ffi^"*/"^ is the rescaled amphtudc of the field and lu'^ ~ ujc/^- In what foUows we assume that 51 = f. 
Furthermore, for simphcity we assume that e is a parameter of the system equal to the initial eccentricity of the orbit 
in the Keplerian problem (wc = and F = Q). 

In this paper, we consider the high-scaled frequency regime for co- rotating orbits : ^/lok > 1 where the Kepler 
frequency lok is equal to i.e., we consider part of phase space corresponding to < J < 1. 



A. Study of the integrable part of Hamiltonian 

Several key features of the dynamics of Hamiltonian (|^) can be understood by looking at the integrable part. We 
consider the mean value with respect to 6 of the integrable Hamiltonian ^ : 

H^^-^-L + J^'^mW (8) 

where = Yln^-oo ' "-"^ another way for considering iJg is to assume that J"*^ X^tit^o m ^mVm-n cosnO is part 

of the perturbation of Hamiltonian . One of the main features of Hamiltonian Ho is that it does not satisfy the 
standard twist condition for luc 0. Since the Hessian of this Hamiltonian is 

the phase space is divided into two main parts : a positive twist region where d^Ho/dJ^ > for J > (\/2/ti;c||I^||)^^'^, 
and a negative one where d^Ho/d.P < for J < (v^/wcHV^H)^^"^. Both regions are separated by a twistless region 
[cPHo/dJ^ = 0). We notice that the singularity at J = is located inside the negative twist region, that the energy 
in the negative twist region is negative, and that the positive twist and twistless regions do not exist in the absence of 

— 1/3 

magnetic field. Furthermore, the size of the negative twist region decreases like ujc as one increases the magnetic 
field uJc, i.e., it shrinks to the singularity of the Hamiltonian (J = 0). 

The phase space of Hq, as well as the one of Hamiltonian (Q), is foliated by invariant tori, the main difference being 
that the invariant tori for Hq are flat in these coordinates. We consider a motion with frequency w, i.e. such that the 
dynamics is 9{t) ~ uot + 0(0) mod 2ti and J{t) = J(0). The associated invariant torus is located at such that the 

frequency dX J = Juj is equal to uj. The equation determining J^^ is : 

2 

There are two real positive solutions of this equation : 



'^ni^ii'^^-^^^ + i = o. (9) 
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The condition of existence of an invariant torus with frequency oj for is lo > V2lJc\\V\\. There are two invariant 
tori with frequency u: one located at in the negative twist region is a continuation of the torus with frequency uj 
in the absence of magnetic field since limj^^^o — the other torus located at in the positive twist region 

is created far from the singularity J = as soon as the field is non-zero. Figure [I] depicts the position of these tori. 
We notice that as soon as Wc 7^ there is creation of a set of invariant tori far from J = in the positive twist region. 
If we increase Wc, the position of the positive twist torus decreases whereas the position of the negative twist torus 
increases. Both tori collide at cUc — uj/y/2\\V\\ to a twistless invariant torus (of the same frequency lj) located at 
{uj/2)~^/^. The approximate location of this twistless torus as a function of uJc (for fixed parameter e) is plotted in 



Fig. N. As one increases uJc, a large portion of the invariant tori with uj G [0, 1] disappears. 
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Remark : First Order Delaunay normalization |lq] 
Averaging Hamiltonian ^ over 6 gives 

{H)e = -^,-L + J'^lWr - I J^Fcos V'. 

Using the expansion of to the second order of e, ~ 1 + and the fact that the previous Hamiltonian 

does not depend on 9 (thus J is constant), the Hamiltonian reduces to 

JC^^L+^e^J^- ^eJ^F cos V, 
16 2 

which is the Hamiltonian studied in Ref. n& to find bifurcation of equilibrium points. 



B. Primary main resonances and Chirikov resonance overlap 

The positive and negative twist regions have their own sets of primary resonances given by Hamiltonian (|^). The 
approximate locations of these primary resonances m:l, denoted J^, are obtained by the condition m9 + ^ w 0. There 
are two such resonances located approximately at 

1/3 



The resonance located at J~ is the continuation of the resonance in the case = 0. The condition of existence of 
these resonances is m\/2tL!c||^|| < 1- Similarly to collisions of invariant tori, collisions of periodic orbits occur when 
increasing tOc- Figure ^ depicts the different domains of existence of real Jm for first primary resonances (m = 1, . . . ,5) 
in the plane of parameters e-ujc- We notice that the most relevant parameter in this problem is the magnetic field. 
The variations of the dynamics induced by the eccentricity e are smooth. 

Between two consecutive resonances m:l and m+l:l (if they exist), regular and chaotic motions occur for a typical 
value of the field F. In order to estimate the chaos threshold between these resonances, i.e. the value for which there 
is no longer any rotational invariant torus acting like a barrier in phase space, Chirikov overlap criterion provides 
an upper bound but usually quite far from realistic values (obtained by numerical integration for instance). For 
quantitatively accurate thresholds, a modified criterion is applied : the 2/3-rule criterion. 

In order to compute the resonance overlap value of the field F between m:l and m+l:l primary resonances, we 
follow the procedure described in Ref. jl^. First, we change the frame to a rotating one at the phase velocity of 
resonance m:l. We apply the following canonical change of variables : {A' , <f') = {N~^A, Ncp) where A = (J, L) and 

= {9, ip) and 



N = 



m 1 
1 



and N denotes the transposed matrix of N. Hamiltonian (|^) is mapped into 

H - -^^^^ - J' "^m^J'^ ^n'Vn'-n cos[n(0' - i,')lm\ 

n.7i' 

+Fm^j'^^Kcos[(n6'' + (m - n)tp')/m]. 

71 

By averaging over the fast variable ip^ the Hamiltonian becomes 

H = -J' + ^^mV^llFf + FmV„ J'2 cos9' . 

2m'' J 8 

The resulting Hamiltonian is integrable and one can compute the width of the resonance m: 1 following Ref. We 
expand the previous Hamiltonian around = Jm/'m and keep only the quadratic part in AJ' — J' — and the 
constant term in the action A J' proportional to cos9' : 

= (l - l^llWrj'^) AJ''+Fj'^V„.cos9'. 
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We rescale energy by a factor — 3rn^J„j'*(l ~ ujIWVW^ : 

1 



H = -A J'- 



F V 



3,„2(1_ 1^2||^||2J6J 



COS 6*'. 



Depending on the positive or negative twist region, this rescaling is positive or negative respectively. 
The width of the resonance m:l in the variables AJ = mAJ' is given by : 



A. 



3 . FVmie) 



-1/2 



since = — 1 + J^/m. The 2/3-rule criterion for the critical threshold for the overlap between resonance 

m:l and m+l:l is reached when the distance between two neighboring primary resonances is equal (up to a factor 
2/3) to the sum of the half- widths of these resonances : 



The critical threshold is given by : 



2 1 

-^\Jm+l — Jm\ = -^{Am + A,„+i). 



3 1 Jrn^^m 



-1/2 



+1 



m-\-l 



-1/2 



(11) 



where stands for either or given by Eq. Therefore we obtain two critical couplings : one in the 

positive twist region, F+(e,ajc), and one in the negative twist region, F~(e,ujc)- We notice that for Wc = 0, since 
J~ = m^/^, we obtain the formula of Ref. for the chaos threshold in the CP problem. For Uc = l/[-\/2\\V\\(rn + l)], 
the threshold Fm{e,ujc) vanishes. This case corresponds to the collision of the resonances m+l:l (the positive and 
negative twist ones). Therefore, the Chirikov criterion is valid only for uj < l/[^/2\\V\\{m. + 1)], i.e. when the two 
primary resonances m:l and m-|-l:l exist in the system. 



We use this criterion in order to study the stability in the different regions of phase space as a function of the 
magnetic field ujc (for small values of the field tOc) and the eccentricity of the initial orbit e. However, since this 
criterion is purely empirical, we use another method to validate or refine the results : we use the renormalization 
method which has proved to be a very powerful and accurate method for determining chaos thresholds in this type of 
models js], |l8[ |l^, 0. We compare the results given by both methods in the region where the criterion applies and 
we use the renormalization map to compute chaos thresholds when Eq. ([ll| ) does not apply (when one of the main 
primary resonance m:l has disappeared). 



C. Renormalization method 



The Chirikov resonance overlap criterion gives us a value for the chaos threshold between two neighboring primary 
resonances. This value aims at approximating the value of the field F for which there is no longer any barrier in phase 
space, and for which large-scale diffusion of trajectories occur between these resonances. The renormalization method 
gives a more local and more accurate information. This method computes the threshold of break-up of an invariant 
torus with a given frequency u. Then by varying w, we obtain the global information on the transition to chaos. 



1. Expansion of the Hamiltonian around a torus 



In order to apply the renormalization transformation as defined in Refs. for a given torus with frequency 

UJ, we expand Hamiltonian (^) in Taylor series in the action J around J ~ J^^. 



-— - Vc_i-)fc+i "-'^-^ A/fc 



k + 1 

fc=0 

where A J ~ J — J^^. The meanvalue of the quadratic term in A J of Hq is equal to ^^^^^jJ^^AJ^. We rescale the 
action variables such that this quadratic term is equal to 1/2, i.e. we replace 77(A J, L, 6*, -0) by Aff (A J/A, L/A, 6*, -0) 
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with A = 5{ujJ^ — 2)/J^. We notice that this rescahng can be done except for two cases J^^ = (a;/2)~-'-/^ which 
is the twistless case, and for J„ = which is the singularity. Therefore, as it is defined in Refs. |5[ |2^, the 
renormahzation cannot be apphed in the twistless region. 
The Hamiltonian (0) becomes : 



+ 00 

E 

fc>2 fc=On>0 fc=On=-oo 



where 



H2,0 



1 

2' 

9(2-L.j3)2' 
J^(11-C.J3) 

108(2- wj3)3' 



3 



m— — oo 



-f^3,n — 



2J3(C^J3-1) E™r-ooK^^^ 



m '171— n 



Hi, 



9(2-C.j3)2 ||l/||2 



m m — n 



54(2-c.j3)3 ii^ii 
2 - J3 



3(2-c.J3)' 

for n G Z. The resulting Hamiltonian is of the form : 

HF{A,ip)=LO-A + V{n-A,ip), (13) 

where a; = (w, —1) and $1 = (1, 0), and the coordinates are A = (A J, L) and (f = {6, For small F, the Kolmogorov- 
Arnold-Moser (KAM) theorem states that if ui satisfies a Diophantine condition, the invariant torus with frequency 
uj will persist. In fact, the picture that emerges from numerical simulations is that if F (in absolute value) is smaller 
than some critical threshold i^c(i^) then Hamiltonian (^ [or equivalently (p^)] has an invariant torus with frequency 
u!. If F is larger than this critical threshold, the torus is broken. In order to compute numerically the critical function 
Lu Fciuj), we use the renormahzation method described briefly below. 



2. Renormahzation method 



Without restriction (up to a rescaling of time), we assume that lo £ [0,1]. The renormahzation relies upon the 
continued fraction expansion of the frequency lu of the torus: 

UJ = — = [ao,ai, . . .]. 
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This transformation will act within a space of Hamiltonians H of the form 

H{A,ip)=Lj ■ A + V{n- A,if), (14) 

where Jl — (l,a) is a vector not parallel to the frequency vector us — (tj,— 1). We assume that Hamiltonian (|lj) 
satisfies a non-degeneracy condition : If we expand V into 



E 

fc>0 



(21 

we assume that Vg is non-zero. This restriction means that we cannot explore the twistless region by the present 
renormalization transformation. 

The transformation contains essentially two parts : a rescaling and an elimination procedure jj] . 

(1) Rescaling : The first part of the transformation is composed by a shift of the resonances, a rescaling of time 
and a rescaling of the actions. It acts on a Hamiltonian H ds H' — H oT (see Ref. |^ for details) : 

H\A, v) = \u-^H (-\NaA, ~N-^^ , (15) 



with 

\^2uj-\a + afV^^' , (16) 

and 

Na^[\l), (17) 

and a is the integer part oi lo^^ (the first entry in the continued fraction expansion). This change of coordinates 
is a generalized (far from identity) canonical transformation and the rescaling A is chosen to ensure a normalization 
condition (the quadratic term in the actions has a mean value equal to 1/2). For H given by Eq. (^^, this expression 
becomes 

H'{A, if)=uj' ■ A + Y^ y^t*--) {W ■ Afe'''-'^, (18) 

where 

= [uj', -1) with J = uj^^ - a, (19) 
rj' = (l,a) with a' = l/(a + a), (20) 

V^^W = rkV[';>^ with rk = i-l)'2^-'^Lu'-\a + a)^-'^ (v^^^^ . (21) 

We notice that the frequency of the torus is changed according to the Gauss map : 

u; u;' = uj-'^ - [uj-^] , (22) 

where denotes the integer part of ui^^. Expressed in terms of the continued fraction expansion, it corresponds 

to a shift to the left of the entries 

UJ = [ao,ai,a2, . . .] i-> w' = [ai, 02,03, . . .]. 

(2) Elimination : The second step is a (connected to identity) canonical transformation U that eliminates the 
non-resonant modes of the perturbation in H'. Following Ref. |22], we consider the set C Z'^ of non-resonant 
modes as the set of integer vectors u = (i'i,V2) such that \v2\ > The canonical transformation hi is such that 
H" = H' olA does not have any non-resonant mode, i.e. it is defined by the following equation : 



I-{H' oU) = 0, 



(23) 
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where I is the projection operator on the non-resonant modes; it acts on a Hamiltonian (|lj) as 



kXi 



We solve Eq. (23) by a Newton method foUowing Refs. p3[ . 

Thus the renormahzation acts on H for a torus with frequency w as H" = TZ{H) = H o T o U for a torus with 
frequency uj' . 

The critical thresholds are obtained by iterating the renormahzation transformation TZ. The main conjecture of the 
renormahzation approach is that if the torus exists for a given Hamiltonian H, the iterates TV^H of the renormahzation 
map acting on H converge to some integrable Hamiltonian Hq. This conjecture is supported by analytical results 
in the perturbative regime , and by numerical results |l^, 20 . For a one-parameter family of Hamiltonians 



{Hp}, the critical amplitude of the perturbation Fc{uj) is determined by the following conditions 

7^"i^F ^ Ho{A) = a; • A + i(0 • A)2 for F < F^{uj), (24) 

n— *oo 2 

n"HF oo forF>Fc(w). (25) 

n — *oo 

The critical threshold Fc{uj) is determined by a bisection procedure. 

In order to obtain the value Fm of chaos threshold between resonances m:l ans m-|-l:l, we vary uj between 1/ (m -|- 1) 
and 1/m. The critical threshold is given by Fm — ^^^uje[i/{m+i),i/m] Fc{w). 



III. NUMERICAL COMPUTATION OF CHAOS THRESHOLDS 



A. Chaos thresholds in the negative twist region 

Figure ^ represents a typical plot of the critical function uj i— > F~ (lo) in the negative twist region obtained by the 
renormahzation method for e — 1 and Uc = 0.3. This function vanishes at (at least) all rational values of the frequency 
u! since all tori with rational frequency are broken as soon as the field is turned on. The condition of existence of a 
torus with frequency u; obtained from the integrable case (iTp is w > \/2wc||T/||, which is in that case w > 0.67. We 
expect the tori with frequency belonging to [0,0.67] to be broken by collision with invariant tori in the positive twist 
region before this value of the field ujc, as it is the case in Hamiltonian (^. 

From Fig. H, we deduce that for F > Fc = 0.023, no invariant tori are left in this region of phase space. The 
frequency of the last invariant torus is equal to (7 -I- l)/(27 -|- 1) w 0.7236 where 7 = (\/5 — l)/2. This value of the 
frequency of the last invariant torus varies with the parameters e and ujc (for ojc — 0, see Fig. 2 of Ref. [|9j). The 
main feature of this frequency is that it remains noble as the parameters e and ujc are varied, in the sense that the 
tail of the continued fraction expansion of this frequency is a sequence of 1, or cquivalcntly this frequency expresses 
like (a7 -I- b)/{c"f + d) where a, h, c, d are integers such that ad — be = ±1. 

We compute chaos thresholds between two successive primary resonances by two methods : the 2/3-rule criterion 
and the renormahzation map. We compute the critical value of the field F for which there is no longer any invariant 
torus between resonances 1:1 and 2:1, located at and respectively, as a function of the magnetic field tUc 
and the parameter e in the negative twist region, i.e. F^ (e, Uc) — niaX;jgjx/2,i] P/T coc)- Figure ^ represents two 
computations : for e = 0.5 and for e = 1. 

For LOc — 0, we obtain the values that have been obtained in Ref. The 2/3-rule criterion gives a very good 

approximation for low values of the field lUc (typically for lUc G [0, 0.15]). For small ujc, we develop the critical function 
F-^^{ujc) given by Eq. (pT|). The corrections to F~ due to the magnetic field are of order since J~ = n^^^ + O(w^). 

We notice that at some value of ujc, the curves F-^(e,U!c) given by the Chirikov criterion decrease sharply to zero. 
We have seen that the approximate condition of existence of the resonance m:l, derived from the integrable case (|^), 
is cuc < l/(V2w||F||). For e = 1 and m = 2, this conditio n is uic < 0.23 (and for e = 0.5 this condition becomes 
LUc < 0.30). This means that the criterion derived in Sec. |II B| is inapplicable for ujc larger than these values {J2 
becomes complex) . One way of extending this criterion for larger values of uJc would be to consider the overlap of 
higher order resonances. From the renormahzation map results, we see that after 2:1 resonance disappear, there is 
an increase of stability which can be explained by the fact that the invariant tori in this region of phase space are no 
longer deformed by this neighboring resonance. We expect this region to be fully broken before 1:1 disappears (i.e. 
all the region between 1:1 and 2:1 primary resonances has disappeared), and this happens at uJc ~ 0.45 for e — I, and 
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at bJc ~ 0.6 for e = 0.5. These results are consistent with the sharp decreases of Fc found by renormalization on Fig. |j. 

Furthermore, we investigate the chaos thresholds for the different primary resonances rrvA by the Chirikov criterion. 
Using Eq. we compute the resonance overlap value of m:l and m+l:l for increasing m = 1, . . . , 10. For a fixed 
eccentricity e = 1, Fig. |^ shows the resonance overlap values as a function of the field lOc 

From this figure it emerges that the dynamics in that region is mainly insensitive to the magnetic field for low values 
of this field. The sharp decreases of the curves are due to the disappearance of one of the primary resonances from 
which the resonance overlap is computed. At these values and for larger we expect stability enhancement due to 
the disappearance of this primary resonance. For Wc > l/(V2'7z||F||), the critical curve F~ sharply decreases to zero 
due to the disappearance of all the region between primary resonances m:l and m+l:l into the twistless region. For 
instance, for m = 3, the critical curve F^{uJc) increases slowly with for lj^ G [0,0.11], then for ujc S [0.11,0.15], we 
expect stabilization enhancement by the magnetic field, and for ujc larger than 0.15, the region of phase space between 
3:1 and 4:1 resonances disappears into the twistless region and F^{ujc) vanishes. 

The region between resonances of low order (to small) seem to be more stable for e = 1. However, this feature 
varies with the parameter e as it has been observed in Ref. jl^. For e S [0, 0.8] and for loc ~ 0, the regions near to:1 
with large to become very stable (and this stabilization is increased by the field for low values of Wc) and the diffusion 
of the trajectories is very limited. Therefore, the orbits with eccentricity close to one are the easiest to diffuse in a 
broad range of phase space. In particular, these orbits can ionize more easily than medium eccentricity ones. This 
feature is expected to hold for small Uc (typically for u)c < 0.05) even if the diffusion is now limited in the negative 
twist region since we expect the twistless region to be very stable and act as a barrier in phase space [^ . Figure ^ 
displays the values of F~{e) for e € [0, 1], and for two values of uJc ■ c^c = O.Ol and uJc = O.l. 

For Lu = 0.1, since a large portion of phase space between primary resonances with to large has disappeared into 
the twistless region, the orbits with medium eccentricity become as easy to ionize as the ones with eccentricity close 
to one in the remaining part of the negative twist region of phase space. For eccentricities close to one, the critical 
threshold is dominated by the chaos threshold between resonances 1:1 and 2:1. 

In summary, the effect of the magnetic field Uc is to stabilize the dynamics in the negative twist region by successively 
breaking up primary resonances. Therefore, between two succesive resonances with to small, the effect of the magnetic 
field is expected to be smooth in the sense that the critical thresholds are only slightly changed (increased) from the 
chaos thresholds in the absence of magnetic field up to some critical value of lOc where collisions of periodic orbits 
occur. Investigating this region by Chirikov resonance overlap yields values which are qualitatively and quantitatively 
accurate with respect to the renormalization results for ujc G [0,0.1]. 

For small values of Uc, the qualitative behavior concerning diffusion of trajectories is expected to be the same as 
the one in the absence of magnetic field even if the diffusion coefficients may be smaller due to the limited region of 
phase space and the stability enhancement due to the magnetic field. Increasing loc makes the orbits with medium 
eccentricity easier to diffuse in a more and more reduced region of phase space. 



B. Chaos thresholds in the positive twist region 

The classical dynamics in the positive twist region of phase space is essential for the stochastic ionization process 
since it contains the region where the action J is large, anf for large Uc this region is the predominant region since 
the negative twist region shrinks to the singularity J — 0. 

Figure]^ show a typical critical function uj i— > F^{uj) obtained by the renormalization method for coc — 0.3 and for 
e — I. This figure is analogous to Fig. ^. 

From Fig. ^ we deduce that for F > Fc — 0.021, no invariant tori are left in this region of phase space. For 
these values of the parameters e and uic, the condition of existence of invariant tori ui > ^/2ujc\\V\\ obtained from the 
integrable case (|^), is w > 0.67. We notice that there is a large chaotic zone for u) S [0.67, 0.8] (the set of w where F+(w) 
is very small). The diffusion of trajectories throughout the J large region is prevented by the small set of invariant 
tori with frequency ui S [0.85, 0.90] for intermediate values of the amplitude of the microwaves Fc € [0.015, 0.02]. The 
frequency of the last invariant torus to survive in this region of phase space is approximately equal to 0.87. The 
frequency of the last invariant torus in the positive twist region fluctuates in a very erratic way as one varies uJc- 
However we observed that this frequency is always between 0.85 and 0.90 (close to the 1:1 resonance). 

The changes induced by the magnetic field are stronger in the positive twist region. We compute the critical 
thresholds between 1:1 and 2:1 resonances in the positive twist region. Figure ^ represents the results obtained by 
renormalization and by the 2/3-rule criterion. In this region, the critical thresholds determined by 2/3-rule criterion 
are well below the ones given by the renormalization for significant values of the field. 

Figure ^ represents the critical threshold for the break-up of invariant tori with frequency 7 — (VS — l)/2. In the 
negative twist region, the golden mean torus is slightly stabilized by the magnetic field. In the positive twist region, 
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the influence of the field is very strong : first, the torus is created then stabihzed to high values of the field, then 
it disappears. This figure shows that in contrast to the negative twist region, the infiuence of the magnetic field in 
the positive twist region is very strong. The curve for the positive twist torus appears to have non-smooth variations 
conversely to the negative twist torus like for instance for uJc ~ 0.065. We notice that for the golden mean tori, there 
is no break-up by collision since the one located in the positive twist region is broken before the expected collision. 

For small values of the field uJc, we expand the threshold given by Chirikov criterion (|l^). Since the resonances 
located at J+ do not exist in the absence of the field lOc, we expect F^{m, e, uJc) to vanish at uic = 0. The corrections 
to F+ due to the magnetic field are obtained using the following expansion 

Thus we have : F+(e, Wc) = Q!m(e)wc^^[l + 0{uj'^/^)\, where am are some constants depending on e. This means that 
the increase of stabilization is very sharp (with infinite slope) for low values of lOc- 

Furthermore, using Eq. ([Tl|), we compute the 2/3-rule resonance overlap value of m:l and m+l:l resonances for 
increasing m = 1, . . . , 10. For a fixed eccentricity e ~ \ , Fig. ^ shows the resonance overlap values as a function of 
the field uJc in the positive twist region. 

What emerges from this figure is that this region of phase space is strongly stabilized by the magnetic field, and 
that the curve m ~ 1 dominates. However, this feature may vary with the parameter e as it has been observed in 



Ref. and in Sec. pi A| for the negative twist region. 

These curves are similar to Fig. ^. As a result, increasing the magnetic field makes intermediate eccentricity orbits 
as easy to ionize (diffusion through a part of phase space where J is large) as the ones with eccentricity close to one. 
This situation is different from the situation without magnetic field. 

The overall effect of the magnetic field in the positive twist region is basically the same as the one in the negative 
twist region, i.e., it breaks up the primary resonances and fills up the remaining region by very stable invariant tori, 
preventing the diffusion of trajectories throughout phase space. 



IV. CONCLUSION 



We have analyzed the classical phase space of the hydrogen atom in crossed magnetic and microwave fields in the 
high frequency regime. Useful information about the dynamics is provided by the analysis of an integrable part of 
the Hamiltonian. Accurate information about chaos threshold is obtained by using the renormalization map and 
the 2/3-rule Chirikov criterion. The global effect of the magnetic field is to stabilize the dynamics and consequently 
reducing the diffusion of trajectories and the stochastic ionization process. 
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FIG. 1: Position of the invariant tori with frequency (\/5 ^ l)/2 (continuous curve) and (5 + V5)/10 (dashed curve) as a 
function of tJc for the integrable Hamiltonian (^) for e = 1. The strong continuous curve is the location of the twistless region. 
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FIG. 2: Existence of primary resonances (m = 1, . . . , 5) in tlie plane of parameters e-Wc- The light gray part is the domain of 
existence of only m = 1, and in the black part, all the five first resonances exist. 
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FIG. 3: Critical function (a») in the negative twist region, obtained by the renormalization method, for ujc = 0.3 and for 
e= 1. 
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FIG. 4: Chaos thresholds F-^ (wc) between resonances 1:1 and 2:1 in the negative twist region, obtained by the 2/3-rule 
criterion (dashed curves) and by the renormalization method (continuous curves), as a function of uJc, for e = 0.5 and e = 1. 
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FIG. 5: Resonance overlap values of F~(wc) between resonances m:l and m+l:l for m = 1, . . . , 10 in the negative twist region 
for e = 1. 
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FIG. 6: Resonance overlap values of F^{e) between resonances m:l and m+l:l for m = 1, . . . in the negative twist region for 
(a) iVc = 0.01 and (6) iVc = 0.1. 
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FIG. 7: Critical function in the positive twist region, obtained by the renormalization method, for ujc = 0.3 and for 

e= 1. 
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FIG. 8: Chaos thresholds F^{u!c) between resonances 1:1 and 2:1 in the positive twist region, obtained by the 2/3-rule criterion 
(dashed curves) and by the renormaUzation method (continuous curves), as a function of ujc- The upper curves are obtained 
for e = 1 and the lower ones are for e = 0.5. 
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FIG. 9: Critical threshold Fc("/,lUc) for the break-up of the invariant tori with frequcucy 7 = (V^ - l)/2 for e = 1 in the 
positive twist region (continuous curve) and in the negative twist region (dashed curve). 
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FIG. 10: Resonance overlap values of F^{uic) between resonances m:l and m+l:l for m = 1, . . . , 10 in the positive twist 
region for e = 1. 
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FIG. 11: Resonance overlap values of F^{e) between resonances m:l and m+l:l for m = 1, . . . in the negative twist region 
for (o) ujc = 0.01 and (6) Wc = 0.1. 



